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Previously  we  showed  that  the  thermoelectric  (TE)  performance  of  bulk  n-type  Bi2Te2  7Se03  can 
be  enhanced  by  subjecting  it  to  a  combined  process  of  chemical  or  mechanical  exfoliation  (C/ME) 
followed  by  a  rapid  densification  and  restacking  of  the  exfoliated  layers  via  the  spark-plasma-sintering 
technique  (SPS).  Here,  we  present  a  systematic  micro-Raman  study  of  two-dimensional  flakes  of  n- 
type  Bi2Te27Se03  produced  by  the  C/ME  process,  as  a  function  of  the  flake  thickness.  We  found  Raman 
evidence  for  flakes  with:  (i)  integer  number  of  quintuples  which  exhibited  a  strong  electron-phonon 
coupling,  and  (ii)  non-integer  number  of  quintuples,  or  sub-quintuples  which  exhibited  the  forbidden 
IR  active  mode  due  to  symmetry  lowering.  Detailed  atomic  force  microscopy  was  used  to  confirm  the 
number  of  quintuples  in  all  flakes  examined  in  this  study.  The  restacking  and  densification  of  these 
flakes  by  SPS  promoted  the  formation  of  charged  grain  boundaries,  which  led  to  the  enhanced  TE 
properties  via  the  energy  filtering  process. 


Bulk  pristine  (undoped)  and  doped  Bi2Te3  are  some  of  the  most  efficient  room  temperature  thermoelectric  (TE) 
materials  for  sustainable  power  generation  and  refrigeration  applications1-3.  The  efficiency  of  a  TE  material  is 
determined  by  its  dimensionless  figure  of  merit,  ZT  =  a2crT/Ky  where  a  is  Seebeck  coefficient  or  thermopower, 
a  is  electrical  conductivity,  and  k(  =  k,e  +  kl)  is  the  total  thermal  conductivity,  which  is  comprised  of  electronic 
and  lattice  contributions,  respectively  The  main  challenge  for  improving  ZT  of  any  TE  material  is  the  inherent 
coupling  between  ay  cr,  and  k  that  prevents  increased  a  and  a  with  a  concomitant  decrease  in  k.  Bi2Te3  is  also  of 
tremendous  interest  as  a  topological  insulator  (TI),  which  in  turn  is  promising  for  TE  energy  conversion4-6.  In 
contrast  to  ordinary  materials,  backscattering  of  electrons  due  to  collisions  with  impurities  and  defects  in  the  crys¬ 
tal  lattice  is  completely  suppressed  on  the  surfaces  of  TIs,  giving  rise  to  improved  charge  transfer  and  mobility,  and 
thus  enhanced  a.  The  existence  of  surface  states  in  TIs  arise  from  the  intrinsic  spin-orbit  coupling  that  is  enhanced 
with  increasing  atomic  masses,  which  can  also  increase  the  mass  fluctuation  scattering  effect  to  reduce  kl. 
It  is  thus  observed  that  most  topological  insulators  such  as  Bi2Te3,  Bi2Se3,  exhibit  excellent  TE  properties. 

In  the  1990s,  Hicks  and  Dresselhaus  predicted  that  an  increase  in  the  density  of  states  could  lead  to  enhanced  a7. 
In  addition  to  improving  electronic  transport  in  low  dimensional  materials,  nanostructuring  via  ball-milling  and 
melt-spinning  have  been  effective  in  reducing  k,l  through  increased  phonon  scattering  effects,  resulting  in  an 
enhanced  ZT  in  TE  nanomaterials8,9.  In  recent  years,  chemical/mechanical  exfoliation  (C/ME)  of  layered  materi¬ 
als  has  enabled  the  fabrication  of  two-dimensional  (2D)  nanosheets  that  exhibit  superior  TE  properties  compared 
to  their  bulk  counterparts10-14. 

Bulk  Bi2Te3  exhibits  a  layered  crystal  structure  and  its  conventional  unit  cell  is  comprised  of  three  quintu¬ 
ples.  Each  quintuple  contains  five  atomic  layers  with  a  sequence  of  Te1-Bi-Te2-Bi-Te1  called  quintuple  layers 
(QL),  and  the  weak  van  der  Waals  bond  between  Te^Te1  couples  the  quintuples  together15.  Therefore  during 
the  C/ME  process,  the  Te^Te1  bond  between  the  QLs  can  be  broken,  giving  rise  to  unique  optical  and  elec¬ 
tronic  transport,  compared  to  the  bulk16,17.  Previously  we  employed  the  combined  technique  of  C/ME  followed  by 
spark-plasma-sintering  (SPS),  which  led  to  preferential  scattering  of  electrons  at  charged  grain  boundaries,  and 
significantly  improved  the  TE  compatibility  factor  and  stabilized  the  ZT  peak  at  higher  temperatures  (350-500  K) 
in  n- type  Bi2Te2  7Se0  310.  The  C/ME-SPS  process  resulted  in  two  important  contributions:  (i)  an  increase  in  the 
electrical  conductivity  due  to  an  increase  in  carrier  concentration  despite  the  presence  of  numerous  grain  bound¬ 
aries,  and  (ii)  the  mitigation  of  the  bipolar  effect  via  band  occupancy  optimization  which  led  to  an  upshift  and 
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Figure  1.  X-ray  diffraction  patterns,  (a)  Bulk  n-type  Bi2Te2  7Se0  3  sample  compared  to  the  8-hr  exfoliated  flake 
before  and  after  SPS  with  a  Bi2Te3  quintuple  as  inset;  (b)  (006)  peak  suppressed  in  the  exf 8h  (X100),  but  more 
pronounced  in  the  exf8/z-SPS. 


stabilization  of  the  ZT  peak  over  a  broad  temperature  range  of  -150  K10.  Both  (i)  and  (ii)  implied  the  creation  of 
charged  grain  boundaries  in  Bi2Te2  7Se0  3  due  to  the  C/ME-SPS  process. 

To  understand  the  structural  changes  in  the  Bi2Te2  7Se0  3  brought  about  by  the  C/ME  process,  we  performed 
a  detailed  micro-Raman  study  of  C/ME  n- type  Bi2Te2  7Se0  3  with  varying  layer  thicknesses  down  to  -2  nm  (-2 
quintuples).  To  the  best  of  our  knowledge,  previous  Raman  studies  of  bulk  and  few-layered  Bi2Te3  delineated  the 
Raman-active  modes  in  pristine  Bi2Te3,  and  the  effect  of  Se  doping  in  the  chemically  exfoliated  layers  was  still 
largely  ignored.  Our  transport  measurements  were  performed  on  the  exfn/z-SPS  samples  (exf nh,  where  n  repre¬ 
sents  the  exfoliation  time  in  hours)  and  showed  a  decrease  in  electrical  resistivity  concomitant  with  an  increase  in 
carrier  concentration,  thereby  resulting  in  a  lower  Seebeck  coefficient.  Moreover,  we  found  a  decrease  in  thermal 
conductivity  with  decreasing  layer  thickness.  Micro-Raman  measurements  on  the  C/ME  samples  revealed  the 
co-existence  of  whole  quintuples  with  strong  surface  states  (electron -phonon  coupling)  and  sub -quintuples  that 
exhibited  forbidden  (IR-active)  modes  in  the  Raman  spectra.  The  combination  of  these  two  “phases”  results  in  the 
formation  of  charged  grain  boundaries  upon  SPS-processing  and  thus  enhanced  TE  performance. 

Results  and  Discussion 

Bulk  Bi2Te3  exhibits  a  trigonal  crystal  structure  belonging  to  the  space  group  R3m18,  but  is  more  commonly  rep- 
resented  by  a  hexagonal  crystal  structure  (see  Supplementary  Fig.  SI).  The  highly  anisotropic  hexagonal  unit  cell 
consists  of  three  quintuples  each  consisting  of  five  atoms,  stacked  in  the  order  Te1-Bi-Te2-Bi-Te1  along  the  c-axis, 
with  lattice  constants  a  =  4.38  A  and  c  =  30.36  A1.  Each  quintuple  measures  approximately  1  nm  across  the  five 
atoms.  The  quintuples  are  held  together  by  weak  van  der  Waals  forces  (Te^Te1  bond)  corresponding  to  the  largest 
spacing  d~  0.37  nm16  that  make  them  easily  cleavable.  In  general,  for  Bi2Te3_xSex  the  Se  atoms  preferentially 
replace  Te  at  Te2  sites  first  and  then  randomly  replace  Te  at  the  Te1  sites14,19.  With  Se-doping  at  the  Te2  sites,  a 
breakdown  of  the  Bi2Te2  7Se0  3  quintuples  into  bi-layer  or  tri-layer  sub -quintuples  (e.g.,  Bi-Te1,  Te^Bi-Te2, 
Te^Bi-Se2)  during  the  C/ME  process  is  feasible  because  the  Bi-Te1  bond  strength  is  the  strongest  bond  in  the 
quintuple20.  In  addition,  the  Te2  atom  is  known  to  lie  at  the  inversion  center  of  the  D3/  symmetry21,  and  hence  the 
Se  doping  can  change  the  crystalline  structure  as  well  as  the  lattice  dynamics  in  a  unique  manner. 

Figure  la  shows  the  powder  HR-XRD  patterns  of  the  bulk  compared  to  that  of  the  exf 8h  and  exf8/z-SPS  sam¬ 
ples.  Both  XRD  patterns  are  consistent  with  the  expected  pattern  for  Bi2Te2  7Se0  3  (JCPDS  card  no.  00-050-0954). 
The  (006)  peak  (Fig.  lb)  was  severely  broadened  in  the  data  of  the  exf 8h  sample  (inferred  by  magnifying  it  a  100 
fold)  due  to  nano-structuring  and  bond  cleavages  at  various  locations  along  the  c-direction,  while  it  became 
more  pronounced  after  SPS  treatment,  suggesting  an  improvement  in  the  coherence  length.  The  representative 
TEM  and  AFM  images  (Fig.  2a  and  b)  of  the  exf 8h  samples  show  lateral  sizes  ranging  from  -0.3  fim  to  0.8  ji m 
and  a  height  of  -2  nm  (Fig.  2d).  The  roughness  (R^)22  of  this  representative  AFM  image  is  about  0.1  nm.  The 
average  thickness  of  the  n- type  Bi2Te2  7Se0  3  layers  as  a  function  of  the  exfoliation  time  is  shown  in  Fig.  2c.  With 
an  increase  in  the  exfoliation  time,  the  layer  thickness  t  of  the  Bi2Te2  7Se0  3  decreased  significantly  between  the 
exfoliation  times  of  0  to  3  hr,  down  to  - 1.6  nm  after  8  hr  exfoliation.  Further  exfoliation  time  did  not  significantly 
reduce  the  thickness  of  the  flakes.  The  presence  of  sub -quintuples  was  evidenced  from  the  non-integer  values  of  t 
from  AFM  linescans  of  the  flakes,  since  1  nm  corresponds  to  one  quintuple  as  discussed  above. 

Figure  3  shows  the  thermoelectric  transport  properties  of  the  exfn/z-SPS  samples  compared  to  the  bulk 
sample.  The  electrical  resistivity  ( p=  1/a)  increased  linearly  with  increasing  temperature  (Fig.  3a),  indicating  a 
“metal-like”  or  degenerate  semiconducting  behavior.  A  slight  change  of  slope  was  observed  in  the  bulk  sample  at 
-360  K,  consistent  with  the  upturn  in  a  (Fig.  3c)  due  to  the  bipolar  effect.  With  increasing  exfoliation  time,  the  p 
values  in  the  exfn/z-SPS  samples  were  consistently  lower  than  the  corresponding  value  in  the  bulk  due  to  increase 
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Figure  2.  Representative  TEM  (a)  and  AFM  (b)  images  of  the  exf 8h  sample,  (c)  The  average  thickness  of  the 
exfoliated  n- type  Bi2Te2  7Se0  3  flakes  as  a  function  of  exfoliation  times  (0  h,  3  h,  5  h,  and  8  h).  The  full  width  of  half 
maximum  of  thickness  distribution  (Gaussian  distribution)  were  used  as  error  bars,  t  tends  to  be  infinite  at  Oh 
and  the  dash  line  is  a  guide  to  the  eye.  (d)  The  lateral  size  and  height  of  a  representative  flake  shown  in  (b). 


in  the  carrier  concentration  n  (Fig.  3b  inset).  The  carrier  mobility  p  —  1/nep,  where  e  is  the  electron  charge, 
(Fig.  3b)  also  exhibited  a  reduction  with  increased  exfoliation  times  at  temperatures  below  100  K  where  the  grain 
boundary  and  charged  defect  scattering  dominates.  Puneet  et  al.  attributed  the  increase  in  n  to  selective  filtering 
of  charge  carriers  by  positively  charged  grain  boundaries  resulting  from  Te  Bi  anti-site  and  bismuth  vacancies  (VBi 
defects)10,16. 

The  Seebeck  coefficient  or  thermopower  (a)  of  the  bulk  Bi2Te2  7Se0  3  sample  (Fig.  3c)  exhibited  a  high  negative 
value  (-200  /xV/K)  at  300  K,  while  a  in  the  exfn/z-SPS  samples  reduced  to  130-170  /xV/K  (65-85%  of  the  bulk 
value),  which  is  consistent  with  the  reduction  of  p  due  to  increase  in  n ,  as  a  varies  as  n~213  shown  in  equation  (1). 
Using  a  simple  parabolic  band  approximation,  the  effective  masses  (m*)  of  the  charge  carriers  in  exfn/z-SPS  sam¬ 
ples  were  calculated  from  the  measured  a  and  n  values  at  room  temperature  using  the  relation23,24. 


StvX 


3  eh2 


m*T\  — 
3  n 


(2/3) 


(1) 


where  kB  is  the  Boltzmann  constant,  e  the  electronic  charge,  and  h  the  Planck’s  constant.  The  effective  masses  (m*) 
for  the  bulk  and  exf nh  samples,  estimated  from  a  vs.  n  at  300  K  (Pisarenko  plot,  Fig.  3c  inset)  were  found  to  be 
-0.95  me ,  where  me  is  the  mass  of  the  electron.  In  addition,  the  band  gaps  ( Eg )  of  the  exfoliated  samples  estimated 
using  the  Goldsmid  Sharp  relation25: 


1 -“g  z'CLXmax1max 


(2) 


remained  constant  -0.16  eV  (see  Supplementary  Table  SI)  which  are  consistent  with  the  reported  values1.  This 
indicates  that  the  systematic  reductions  in  a  with  increasing  exfoliation  times  could  be  solely  attributed  to 
increasing  n  (inset,  Fig.  3b)  and  not  from  changes  in  band  gap  or  band  curvature.  With  increasing  temperature, 
|a|  of  the  bulk  sample  increased  till  amax  was  reached,  above  which  the  bipolar  (two  carrier)  effect  was  observed. 
Specifically,  \a\  decreased  as  indicated  by  the  upturn  at  T  =  362  K  in  Fig.  3c,  which  is  typical  of  a  narrow  band 
Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 
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Figure  3.  Temperature  dependence  of  transport  properties,  (a)  The  electrical  resistivity  (p);  (b)  mobility  (/x) 
with  carrier  concentration  (n)  as  inset;  (c)  temperature  dependence  of  the  Seebeck  coefficient  (a)  with  its 
dependence  on  carrier  concentration  shown  in  the  inset;  (d)  total  thermal  conductivity  (k)  with  electronic 
thermal  conductivity  (%)  as  inset  of  the  exf3/z-SPS,  exf5/z-SPS,  exftffi-SPS  samples  as  a  function  of  temperature 
compared  to  the  bulk  sample. 


gap  semiconductor.  As  concluded  by  Puneet  et  al.  the  low  energy  minority  carriers  in  the  exf 8h  samples  were 
selectively  filtered  due  to  the  presence  of  the  positively  charged  defects  on  the  grain  boundary  (interfacial  charged 
defects)  which  upshifted  the  amax  to  above  500  K10.  With  the  reduction  of  the  bipolar  term,  this  upshift  in  amax  led 
to  the  broadening  of  ZT  in  the  exftffi-SPS  samples  as  shown  in  Supplementary  Fig.  S2. 

The  presence  of  grain  boundary  potential  barrier  scattering  in  the  n-type  Bi2Te2  7Se0  3  samples  is  evidenced 
from  the  temperature  dependence  of  cr,  where  a  ~  T_1/2exp(^f  )26  and  the  associated  grain  boundary  potential 

barrier  height  (EB)  was  calculated  from  the  linear  plot  of  ln(crT112)  vs.  1/kT  as  shown  in  Supplementary  Fig.  S3.  As 
expected,  the  exfnh- SPS  samples  exhibited  higher  Es  ~  23-27  meV  arising  from  the  increased  number  of  grain 
boundaries  (grain  size  ~20  /x  m)  compared  to  the  EB~  19meV  of  the  ingot  that  exhibited  relatively  larger  grains. 
For  comparison,  an  EB  ~  60  meV  was  reported  in  PbTe  nanocomposites,  wherein  grain  boundary  potential  barrier 
scattering  is  known  to  be  a  dominant  scattering  mechanism27.  Furthermore,  evidence  for  the  presence  of  charged 
grain  boundaries  in  n-type  Bi2Te2  7Se0  3  samples  was  also  found  using  Kelvin  probe  Force  Microscopy  (KPFM; 
Supplementary  Fig.  S4)  as  detailed  in  the  SI  section. 

The  bulk  sample  exhibited  a  high  magnitude  of  total  thermal  conductivity  (k)  at  ~  17  K  as  shown  in  Fig.  3d, 
which  decreased  subsequently  in  the  exfnh- SPS  samples  due  to  exfoliation-induced  disorder.  A  well-defined 
peak  is  expected  at  lower  temperatures  in  the  bulk  sample  (indicating  good  crystal  quality),  before  k,  decreased 
with  increasing  temperature  due  to  anharmonic  phonon-phonon  or  Umklapp  scattering  effects.  At  -360  K,  k, 
increased  gradually  with  increasing  temperature,  due  to  contribution  from  the  bipolar  effect.  The  inset  showed 
the  systematic  increase  in  ke  ( =L0aT )  (where  L0  is  the  Lorenz  number,  L0  =  1.66  x  10_8V2/K2  for  the  nanos- 
tructured  samples28)  in  the  exfnh- SPS  samples  compared  to  the  bulk,  indicating  that  kt—  ke  was  also  reduced 
due  to  increased  phonon  scattering  at  the  grain  boundaries.  The  bipolar  contribution  however  was  absent  and 
possibly  shifted  to  higher  temperatures  in  the  exfoliated  samples.  Taken  together,  the  plots  in  Fig.  3  show  that 
the  C/ME-SPS  treatment  of  n-type  resulted  in:  i)  an  increase  (decrease)  in  carrier  concentration  (resistivity),  ii) 
mitigation  of  the  bipolar  effect  in  thermopower,  and  iii)  a  simultaneous  reduction  in  the  thermal  conductivity, 
that  led  to  the  broadening  of  the  ZT  peak  over  a  wider  range  of  temperature  -100  K  (see  Supplementary  Fig.  S2). 

To  gain  more  understanding  into  the  structural  properties  underpinning  the  unique  TE  performance  of  the 
exfoliated  samples,  we  used  micro-Raman  spectroscopy,  which  is  an  ideal  technique  for  studying  the  changes 
in  the  vibrational  (and  hence  TE)  properties  caused  by  the  C/ME  process  in  n- type  Bi2Te2  7Se0  3.  Micro-Raman 
spectra  were  collected  from  several  different  spots  on  the  bulk,  exf3fi,  exf 5h  and  exf 8h  Bi2Te2  7Se0  3  samples.  Upon 
densifying  these  samples  using  SPS,  the  SPS-compacted  samples  exhibited  a  high  background  making  them 
unsuitable  for  Raman  studies.  Nevertheless,  significant  structural  changes  in  the  exfoliated  layers  due  to  the  C / 
ME  process  were  evident  in  our  micro-Raman  spectra,  as  discussed  below. 
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Figure  4.  Vibration  modes  and  Raman  measurement,  (a)  Schematic  diagram  of  the  four  Raman-active  modes 
and  two  IR-active  modes  of  the  Bi2Te3  quintuple.  The  purple,  blue  and  green  colors  represent  the  Bi,  Te2  and 
Te1  atoms  respectively.  The  vertical  arrows  represent  out  of  plane  vibration.  The  horizontal  arrows  represent 
in  plane  vibration  that  parallel  to  the  page.  The  — ,  +  signs  represent  in  plane  vibrations  that  are  perpendicular 
to  the  page.  (b,c)  Micro-Raman  spectra  of  n-type  Bi2Te2  7Se0  3  without  and  with  the  presence  of  the  Alu2  mode 
as  function  of  different  exfoliation  times  compared  to  the  bulk  sample.  The  open  symbols  represent  the  raw 
data  and  the  solid  lines  through  the  raw  data  are  the  fits.  The  individual  Lorentzian  fits  are  shown  below  each 
spectrum.  The  red  bars  indicate  the  position  of  the  optical  phonon  peaks  of  bulk  n-type  Bi2Te2  7Se0  3.  The  gray 
(green)  bar  in  (c)  indicate  the  peak  position  of  the  IR-active  Alul  (Alu2)  mode  in  bulk  Bi2Te3. 


Bulk  n-type  Bi2Te2  7Se0  3  is  known  to  exhibit  four  signature  Raman-active  optical  phonons  as  shown  in  Fig.  4. 
Of  these  four  modes,  the  lowest  frequency  E^1  mode  was  difficult  to  resolve  owing  to  the  increasing  spectral 
background  below  40  cm-1.  The  other  three  modes21  are  centered  at  around  62.4  cm-1  (A^1),  102.6  cm-1  (Eg2) 
and  136.5  cm-1  (Alg2)  as  represented  by  the  red  vertical  bars  in  the  spectra  in  Fig.  4b  and  c.  Two  types  of  Raman 
spectra  were  observed  at  various  spots  from  all  samples  irrespective  of  the  exfoliation  time:  those  that  exhibited 
the  three  modes  mentioned  above  (Fig.  4b),  and  those  that  exhibited  an  additional  peak  at  -122  cm-1  (Fig.  4c). 
The  additional  peak  is  an  IR-active  mode  (Alu2)  that  has  been  reported  previously16,29-32,  in  nanoscale  pristine 
Bi2Te3  (mode  at  around  1 16  cm-1  indicated  by  the  green  vertical  bar  in  Fig.  4c).  We  attribute  the  presence  of  the 
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Alu2  mode  to  symmetry  breaking,  possibly  arising  from  the  disorder  induced  by  exfoliation  and/or  Se-dopant  at 
the  Te  sites  (Te1  and  Te2),  where  Te2  is  the  inversion  center  of  the  crystal  symmetry  (see  Supplementary  Fig.  SI)33. 
As  Se  is  lighter  than  Te,  Alu2  mode  was  found  blueshifted  to  -122  cm-1  compared  to  the  corresponding  peak 
frequency  in  bulk  Bi2Te3  (Fig.  4c),  confirming  the  presence  of  Se.  Several  dozen  spots  were  scanned,  and  Alu2 
mode  was  observed  in  the  Raman  spectra  at  roughly  half  of  these  spots,  as  can  be  seen  in  the  2D  Raman  intensity 
maps  in  supplementary  Fig.  S5.  The  random  occurrence  of  the  Alu2  mode  suggests  variations  in  the  structure  of 
the  exfoliated  sheets,  possibly  from  cleaving  of  the  layers  into  sub -quintuples,  as  discussed  further  below.  The 
other  low  wave  number  IR- active  mode  (Alul)  at  -94  cm-1 34  (mode  indicated  by  the  gray  vertical  bar  in  Fig.  4c) 
in  Bi2Te3  could  not  be  discerned  as  this  peak’s  frequency  is  close  to  that  of  the  E^2  mode.  Moreover,  the  Alul  mode 
frequency  may  have  blueshifted  and  as  a  result  could  be  masked  by  the  Eg2  mode. 

We  next  discuss  the  frequency  and  linewidth  dependence  of  the  modes  shown  in  Fig.  5  as  a  function  of  inverse 
thickness  of  the  samples.  Figure  5a  and  c  (Fig.  5b  and  d)  show  the  mode  frequency  and  linewidth  dependences, 
respectively,  in  the  absence  (presence)  of  the  Alu2  mode  in  the  Raman  spectra.  In  Fig.  5a  and  b,  the  A^1  stretching 
mode  frequency  exhibited  the  least  dependence  on  1/t  while  the  E^2  pinch  mode  frequency  increased  slightly  with 
increasing  1/t.  Their  linewidth  dependences  were  relatively  weaker  than  that  exhibited  by  A^2  and  Alul  modes. 
Alg2  mode  is  most  sensitive  to  1/t:  its  frequency  blueshifted  by  -3  cm-1  in  Fig.  5a  while  it  redshifted  in  Fig.  5b  and 
exhibited  a  significant  sharpening  from  -21  to  7  cm-1  with  increasing  1/t  (decreasing  thickness).  A  blueshift  of 
the  A^2  peak  frequency  with  decreasing  layer  thickness  was  also  reported  by  Zhao  et  a/.35  for  CVD -grown  pristine 
Bi2Te3  (inset  in  Fig.  5a)  although  it  was  accompanied  by  a  broadening  of  the  peak. 

The  discrepancy  between  our  results  and  those  from  Zhao  et  al.  can  be  explained  by  considering  the  effect  of 
Se.  The  A^2  phonon  mode  in  Bi2Te3  exhibits  strong  electron  phonon  coupling  (EPC),  which  becomes  stronger 
with  doping36.  The  strong  EPC  is  also  responsible  for  the  formation  of  a  Kohn  anomaly  at  the  Brillouin  zone 
center  and  the  observation  of  Dirac  fermions  in  the  topological  surface  states37.  For  a  phonon  with  strong  EPC, 
the  anharmonic  contribution  to  phonon  decay  is  dwarfed  by  decay  into  electron  hole  pairs.  However,  in  doped 
systems  where  the  Fermi  level  is  greater  than  the  phonon  energy,  Pauli  blocking  reduces  the  number  of  electron 
states  for  the  phonon  to  decay  into,  resulting  in  a  longer  phonon  lifetime.  This  is  manifested  as  a  sharpening  of  the 
Raman  peak.  In  addition,  doping-induced  change  in  the  Fermi  surface  moves  the  Kohn  anomaly  away  from  the 
center  of  the  Brillouin  zone,  where  Raman  active  phonons  are  probed,  and  consequently  causes  a  stiffening  of  the 
phonon  mode.  Indeed  such  observations  have  been  made  in  doped  graphene  and  metallic  carbon  nanotubes38-40, 
where  the  E^2  mode  (G  peak)  blueshifts  and  sharpens  upon  both  hole  and  electron  doping.  While  the  observations 
of  strong  EPC  and  Kohn  anomaly  in  rc-type  Bi2Te3  have  been  made  only  at  low  temperatures  (<20  K),  it  is  possible 
that  the  C/ME  process  breaks  apart  the  Bi2Te2  7Se0  3  layers  into  Se-doped  quintuples  where  Se  is  substituted  at 
the  Te2  sites  in  the  Bi2Te3  quintuple,  maintaining  the  crystal  symmetry  In  that  case  one  could  expect  an  increase 
in  surface  states  especially  with  decreasing  layer  thickness,  consistent  with  our  observation  of  blueshifted  peak 
frequencies  (Fig.  5a)  and  decreasing  linewidths  of  the  Alg  mode  (Fig.  5c)  with  its  strong  EPC.  This  result  also 
corroborates  the  increase  in  electron  densities  and  formation  of  charged  boundaries  surmised  from  the  transport 
measurements  shown  in  Fig.  3. 

Conversely,  in  the  spectra  where  the  Alw2  mode  was  observed  (Fig.  4c),  the  Alg  peak  exhibited  a  redshift  (by 
-4  cm-1)  and  broadening  (by  -5  cm-1),  opposite  to  the  trends  seen  in  the  spectra  in  which  the  Alw2  mode  was 
absent.  As  mentioned  above,  the  appearance  of  the  Alu2  mode  is  attributed  to  breaking  of  crystal  symmetry  from 
the  C/ME  process.  In  this  case  it  is  possible  that  the  C/ME  process  caused  individual  quintuples  to  fragment  into 
sub -quintuples  and  the  appearance  of  the  Alu2  mode  is  induced  by  the  disorder  caused  by  this  fragmentation. 
Due  to  spot-to-spot  variations  we  did  not  observe  any  dependence  of  the  Alu2  mode  intensity  on  1/t.  The  Raman 
modes  did  exhibit  broadening  (by  -2-4  cm-1)  with  increasing  1/t  in  the  spectra  when  the  Alu2  mode  was  pres¬ 
ent  in  the  spectra  (Fig.  5d).  A  decrease  in  the  sample  thickness  typically  leads  to  broadening  of  peaks,  and  such 
a  broadening  was  indeed  observed  in  the  Alg  peaks  with  decreasing  Bi2Te3  layer  thickness  by  Zhao  et  al.35.  The 
observation  of  disorder- induced  IR  active  modes  supports  our  thermal  measurements  on  the  exfn/z-SPS  samples, 
where  k  reduced  compared  to  the  bulk  value  with  increasing  exfoliation  time. 

Based  on  the  Raman  and  transport  measurements,  a  new  picture  of  the  structural  changes  that  occur  in  the  C  / 
ME  processed  Bi2Te2  7Se0  3  crystal  emerges:  i)  Se  dopants  preferably  substitute  for  Te2  sites  followed  by  Te1  sites,  ii) 
during  the  chemical  exfoliation  process,  the  bond  cleavages  at  various  locations  as  depicted  schematically  in  Fig.  6 
(viz.,  Bi-Se  or  Bi-Te1  bond)  to  form  different  sub -quintuples  such  as  Te^Bi-Se2,  Te^Bi-Te2,  Bi-Te1  and  Bi-Se1,  and 
iii)  the  co-existence  of  quintuples  and  sub -quintuples  in  the  exfoliated  samples  (as  evidence  directly  by  AFM  and 
indirectly  by  micro-Raman  studies  described  in  Figs  2  and  4),  promotes  the  formation  of  charged  grain  bounda¬ 
ries  during  SPS  treatment.  The  charged  grain  boundaries  leading  to  the  transport  properties  discussed  in  Fig.  3. 

Conclusion 

In  this  study  we  presented  systematic  TE  properties  of  few  layered  n-type  Bi2Te2  7Se0  3  produced  by  C/ME-SPS 
process.  The  TE  measurements  showed  an  increase  (decrease)  in  carrier  concentration  (resistivity  and  thermal 
conductivity),  accompanied  by  a  mitigation  of  bipolar  effect  in  a  in  the  exfn/z-SPS  samples.  Systematic  AFM 
and  micro-Raman  studies  of  C/ME  processed  samples  provided  evidence  for  the  co-existence  of  quintuples 
and  sub -quintuples,  which  upon  SPS  process  promote  the  formation  of  charged  grain  boundaries.  Specifically, 
micro -Raman  analysis  revealed  two  types  of  spectra,  which  exhibited  different  frequency  and  linewidth  trends  as 
a  function  of  layer  thickness.  The  first  set  of  spectra  provided  evidence  for  scattering  from  Se-doped  quintuples, 
where  the  increased  EPC  leads  to  stiffening  and  sharpening  of  the  Alg  phonon.  In  the  second  set  of  spectra  the 
disorder-induced  IR-active  mode  Alw2  was  evident,  which  we  attribute  to  the  formation  of  sub -quintuples  caused 
by  the  C/ME  process. 
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Figure  5.  Frequency  and  FWHM  dependence  of  lit.  (a)  Dependence  of  the  A^1,  Eg2,  Alg2  mode  frequencies  as  a 
function  of  lit  when  the  Au  mode  is  absent.  Inset  compares  the  Alg2  mode  frequency  (+)  from  Zhao  et  al35  with 
current  work,  (b)  Dependence  of  the  A^1,  Egi  Alg  and  Alu2  mode  frequencies  as  a  function  of  1  It.  The  dashed 
horizontal  lines  in  (a)  and  (b)  represent  the  peak  frequencies  of  Raman  modes  in  bulk  n- type  Bi2Te2  7Se0  3. 
Change  in  linewidths  (FWHM)  of  the  A^1,  Eg,  Alg  modes  as  a  function  of  lit  (c)  when  the  Alu2  mode  is  absent, 
(d)  when  the  Alu2  mode  is  present. 


Methods 

Sample  preparation.  Pieces  of  n- type  Bi2Te2  7Se03  ingot  (Marlow  Industries,  USA)  were  dispersed  in 
N-methyl-2-pyrrolidinone  (NMP)  with  a  ratio  of  lOg/L  and  sonicated  using  1/8-inch  tip  sonicator  (Branson  250) 
at  20  W  for  0, 3,  5  and  8  hrs  to  obtain  Bi2Te2  7Se0  3  flakes.  Subsequently,  the  supernatant  solution  was  centrifuged  at 
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Figure  6.  Schematic  depiction  of  sub -quintuples  formed  as  a  result  of  the  C/ME  process. 


4000  rpm  for  2  hrs  and  the  resulting  powder  was  washed  several  times  using  deionized  water  to  remove  residual 
NMP  and  then  oven  dried  at  -100  °C.  Next,  the  exf nh  flakes  were  compacted  using  spark  plasma  sintering  (SPS, 
Dr.  Sinter  LabH-515S  system)  by  loading  ~2-3  g  of  exf nh  into  graphite  dies  and  sintered  at  500  °C  for  5  minutes  at 
an  applied  pressure  of  30  MPa  under  a  dynamic  vacuum.  The  resulting  SPS  pellets  were  12.5  mm  in  diameter  and 
2-3  mm  in  thickness,  and  a  density  of -98-99%  of  the  theoretical  density  We  refer  to  the  SPS  densified  exf nh  and 
commercial  n- type  Bi2Te3  ingot  samples  as  exfn/z-SPS  and  bulk,  respectively. 

Characterization.  The  thicknesses  of  the  exf nh  flakes  were  measured  using  non-contact  mode  atomic  force 
microscopy  (AFM,  Model:  AIST-NT  Smart  SPM,  Micromasch  cantilevers  HQ:  NSC14/A1  BS-50).  The  average 
thickness  was  statistically  calculated  from  the  AFM  height  measurements  conducted  on  200-400  flakes  for  each 
exf nh  sample.  In  addition,  two-pass  Kelvin  probe  force  microscopy  (KPFM,  Micromasch  conductive  AFM  probes 
HQ:NSC14/Cr-Au,  scan  rate:  1.0  Hz)  was  used  for  imaging  and  measuring  the  contact  potential  difference  between 
the  AFM  tip  and  the  sample  at  the  charged  grain  boundaries.  For  each  scan  line,  during  the  two-pass  KPFM  meas¬ 
urement,  the  height  profile  was  recorded  as  AFM  topographic  image  and  followed  by  lifting  the  probe  by  30  nm 
above  the  surface  to  measure  the  potential  offset.  The  AIST-NT  image  analysis  and  processing  software  (Version 
3.2.14)  was  used  for  AFM  topographic  and  KPFM  image  analysis.  The  microstructural  and  chemical  analyses  were 
performed  using  the  conventional  transmission  electron  microscopy  (TEM,  Hitachi  H7500)  and  high  resolution 
X-ray  diffraction  (HR-XRD,  RIGAKU  Ultima  IV  diffractometer,  Cu  K a  radiation,  A  =  1.5406  A).  Micro-Raman 
spectroscopy  of  exf nh  samples  was  performed  using  a  633  nm  excitation  in  a  Renishaw  Raman  microscope  equipped 
with  a  lOOx  objective  lens  (600  nm  spot  size).  A  reduced  laser  power  (-100  /xW)  was  used  to  prevent  the  inadvertent 
overheating  of  the  exf nh  flakes  during  the  collection  of  their  Raman  spectra.  Similar  micro-Raman  measurements 
were  performed  on  the  exfn/i-SPS  samples  but  the  shiny  surfaces  of  these  samples  did  not  yield  a  good  Raman  signal. 

All  transport  measurements  were  performed  in  the  perpendicular  direction  to  the  SPS  pressing  direction.  The 
temperature  dependent  (15  to  300  K)  resistivity  and  thermopower  were  measured  quasi- simultaneously  using  a 
4-probe  measurement  technique  which  is  described  in  detail  elsewhere41.  The  commercial  ZEM  (ULVAC-RIKO, 
ZEM-2)  was  used  to  measure  resistivity  and  thermopower  from  300-500  K  under  partial  He-atmosphere.  The 
thermal  conductivity  (kt)  was  measured  from  20-320  K,  using  a  standard  steady-state  technique  on  a  custom 
designed  measurement  system42.  The  high  temperature  thermal  conductivity  was  calculated  using  the  relation 
K,T  =  CpDd;  where  d  is  the  packing  density  of  the  material,  D  the  thermal  diffusivity  and  CP  (~Cy,  for  solids) 
the  specific  heat  capacity  The  heat  capacity  was  measured  using  a  NETZSCH  DSC  404C,  thermal  diffusivity 
was  measured  using  a  NETZSCH  LFA  457  system  and  the  packing  density  was  measured  by  the  Archimedes' 
principle.  Since  the  laser  flash  measures  the  thermal  diffusivity  along  the  SPS  pressure  direction,  whereas  the  low 
temperature  transport  properties  are  measured  along  the  direction  perpendicular  to  SPS  direction  (in  the  plane 
of  pellet),  several  bars  of  our  samples  were  cut  and  re-stuck  together  after  rotating  the  bars  by  90  degrees  using 
JB  Weld  (a  thermally  conducting  and  electrically  insulating  glue),  in  order  to  measure  all  the  properties  along 
the  same  direction.  Both  the  electronic  and  thermal  transport  measurements  in  the  low  and  high  temperature 
regimes  were  in  good  agreement  over  the  entire  temperature  range  of  15-500  K.  The  carrier  concentration  was 
determined  from  10-300  K  by  Hall  coefficient  measurements  using  the  commercial  Quantum  Design  Physical 
Properties  Measurement  System  (PPMS)  under  a  magnetic  field  sweep  of  5  kOe. 
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